Constraining Recent Oscillations in Quintessence Models with Euclid 
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Euclid is a future space-based mission that will constrain dark energy with unprecedented accu- 
racy. Its photometric component is optimized for Weak Lensing studies, while the spectroscopic 
component is designed for Baryon Acoustic Oscillations (BAO) analysis. We use the Fisher matrix 
formalism to make forecasts on two quintessence dark energy models with a dynamical equation of 
state that leads to late-time oscillations in the expansion rate of the Universe. We find that Weak 
Lensing will place much stronger constraints than the BAO, being able to discriminate between 
oscillating models by measuring the relevant parameters to la precisions of 5 to 20%. The tight 
constraints suggest that Euclid data could identify even quite small late-time oscillations in the 
expansion rate of the Universe. 



I. INTRODUCTION 

The first evidence for an accelerating Universe was un- 
covered in 1998 through the luminosity distance-redshift 
relation obtained for type la Supernova [TJ [2] , and has 
been consistently corroborated by subsequent observa- 
tions [21 H] . This phenomenon rapidly became one of the 
most intriguing problems in modern Cosmology. Since 
then, there have been many solutions proposed and some 
have survived the increasingly tight constraints imposed 
by observational data and are planned to be analysed by 
Euclid]^ data [5]. Among those solutions, vifhich include 
changes to the left-hand side (or the geometry /gravity 
part) of Einstein's field equations, as in modified gravity 
(/(i?)) theories [HITj, or the introduction of extra spa- 
tial dimensions, as in braneworld models [5], one of the 
most popular assumes the existence of an exotic form of 
energy in the Universe, characterized by a negative pres- 
sure, known as dark energy. If it exists, then present- 
day observations suggest it accounts for approximately 
75 per cent of the Universe's total energy density, and its 
pressure-to-density ratio (known as the equation of state, 
w) is nearly constant with time and close to —1 [SHT^ 
(see Ref. [IS] for more details). If future observations 
continue to point towards a value of w near —1, then the 
hypothesis of the dark energy being in the form of a cos- 
mological constant. A, will gain strength, demanding for 
more insight on the theoretical problems raised by such 
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a scenario (see Ref. [T^ for a review). But it is also con- 
ceivable that observational data will become more easily 
reproducible with a value for w close to, but not exactly 
equal to, —1, opening the door for alternative hypothesis, 
such as quintessence theories. 

In quintessence theories, the dark energy is assumed 
to be the result of a scalar field, (p, minimally coupled 
to ordinary matter through gravity. The field evolves 
in a potential V{(j)) with a canonical kinetic term in its 
lagrangian, which results in a dark energy equation of 
state w > —1, not necessarily constant. These models 
have already been extensively studied [TSHIS] and, usu- 
ally, it is assumed that the evolution of the scalar field is 
monotonic in a potential that satisfies the typical infla- 
tionary slow roll conditions (see Ref. [I^] for a review on 
inflation), as given by: 



idvy 



(1) 



and 



V d(t>^ 



« 1. 



(2) 



Under these conditions, the dark energy equation of state 
will always be close to —1, and the motion of the field 
becomes greatly simplified |20j , allowing for an analytical 
expression for w as a function of the scale factor a. The 
motion of the field will also depend on the present-day 
values of w and of the dark energy density, fl^, where 
is the ratio of the dark energy density, p^, to the 
universe's critical density, pc- 
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However, the inflationary slow-roll conditions need not 
necessarily apply to quintessence models. This opens up 
the possibility of complex scalar field dynamics, allow- 
ing for a more dynamical evolution of the dark energy 
equation of state, while still being on average close to 
— 1. Here we will be particularly interested in scalar field 
models where late-time oscillations in the expansion rate 
of the Universe can arise due to non-standard scalar field 
dynamics. In this paper, we will focus on two scalar field 
models for which such behavior is possible. Model I as- 
sumes a quadratic potential, where an extra degree of 
freedom, related to the curvature of the potential at the 
extremum, is introduced on the evolution of the scalar 
field [2TH23] . More specifically, we consider only the case 
where the curvature of the potential is such that it en- 
ables the field to oscillate around the stable extremum 

In the case of Model II, we consider power-law poten- 
tials where the minimum is zero. When the scalar field 
is close to the minimum, oscillations begin, and an ef- 
fectively constant equation of state arises, averaged over 
the oscillation period (see and references therein). 

In particular, we consider the case where the field starts 
evolving monotonically, and only close to today it starts 
oscillating around the minimum, guaranteeing that the 
observed cosmic history is not significantly affected [57]. 

In Sec. |lTj we review the dynamics of quintessence 
models and briefly present those on which we focus our 
attention. Then, in Sec. |III[ we present the main objec- 
tives and properties of the Euclid mission, and discuss 
the Fisher matrix formalism that will be used to fore- 
cast both the Weak Lensing and BAO constraints. In 
Sec. |IV( we determine the fiducial values for the param- 
eters of each quintessence model, using the most recent 
data available from the Supernova Cosmology Project, 
and then use the Fisher matrix formalism to determine 
to what accuracy Euclid will be able to constrain the pa- 
rameters of the considered quintessence models. In Sec. 
[V] we present some final remarks, establishing a compar- 
ison between our results and those obtained for the usual 
dark energy linear parameters Wq — Wa- Throughout the 
paper we shall assume the metric signature [—,+,+,+] 
and natural units oi c — H = 8ttG = 1. 



II. QUINTESSENCE MODELS 

The action for quintessence is given by 



(3) 



where R is the Ricci scalar, £,„ and £0 are, respectively, 
the matter and scalar field lagrangians. The field's la- 
grangian takes the form 



1 



(4) 



where V{(j)) is the field's potential. 

The energy-momentum tensor of the field is deter- 
mined by varying the field's individual action with re- 
spect to the metric elements g'"', yielding 



Jg"^ (9^09^0) + F(0) 



(5) 



The line element for a fiat, homogeneous and isotropic 
Universe is given by the Friedmann-Robertson- Walker 
(FRW) metric: 



ds^ = -dt^ + a^{t) {dx^ + dy^ + dz^) 



(6) 



where a is the scale factor, t is the physical time and 
X, y and z are the spatial comoving coordinates. Using 
the FRW metric, and assuming that the field depends 
very weakly on the spatial coordinates, the quintessence 
energy density, p, and pressure, p, are given by 



-L n 
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02-T/(0), (7) 



where the dot represents a derivative in order to the phys- 
ical time t. 

As for the dark energy equation of state, w, it will be 
given by 



w = — = 



(1/2)02 +^(^)- 



(8) 



This equation shows that, in order to have a dark energy 
equation of state close to —1, the field's evolution has to 
be potential dominated, such that <j>^ << V{(j)). This is 
why, generically, quintessence models are associated to 
scalar fields slowly evolving in the respective potentials. 
This can be ensured, for instance, by the slow-roll condi- 
tions given by equations ([I]) and ([2|. 

The scalar field equation of motion is given by the 
Euler-Lagrange equation 



3i?0 



dV 



= 0. 



(9) 



Through this equation we can see that the field will evolve 
in the potential V{(f>) rolling towards a minimum in the 
quintessence potential, while its motion is damped by the 
presence of the Hubble parameter, H. Considering a flat 
Universe, with the FRW metric, the Hubble parameter, 
as a function of the scale factor a, is given by 



(10) 



where Pj. is the Universe's total energy density. We will 
be considering a Universe consisting of pressureless mat- 
ter and a scalar field playing the role of dark energy. This 
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implies that = p,„ + P<p- As a function of the scale -o.b5 
factor, pm = PrnocT^ , where pmo is the present-day value 
of the matter density; the field's energy density will be 
given by in Eq. Q. In a flat Universe, the present- _ogg 
day value of the matter energy density is determined by 
^mo 1 ^ ^^00 , where and VL^q are the ratios of the f 
present-day values of the matter and dark energy densi- 
ties to the critical density, pc- 



Model I 




The first model we consider is that of a scalar field 
rolling close to a stable non-zero minimum of its poten- 
tial. The field is assumed to evolve in a potential that 
satisfies the slow-roll condition of equation ([T]), while 
the other slow-roll condition is somewhat relaxed since 
{l/V)d?V / d(tP' can be large, namely as a result of the 
curvature at the potential's minimum. This model was 
extensively studied in [211, where a set of solutions for 
the evolution of the field was derived, as well as the re- 
spective equation of state, w{a), as a function of the scale 
factor a. These solutions are applicable to a wide array 
of potentials with non-zero minima, and depend on the 
present-day values of w and fi^, and on the curvature 
of the potential at the minimum, controlled by the pa- 
rameter . The latter is given by equation (19) in [22], 
which we reproduce here: 



= 1 



4 V"{'P.) 
3 V{(j},) ■ 



(11) 



where V"{(j)i,), each prime representing a derivative with 
respect to the scalar field, and are, respectively, 

the curvature and potential values at the minimum of 
the potential, 0^. The latter corresponds to the present- 
day value of the dark energy density, = p^o. 
For this model we have V"{(j)i,) > 0, which means that 
— c» < < 1. There are three cases of interest: 
K"^ = I, < < 1 and < 0. The hmiting case of 
— > 1, therefore with zero curvature, was considered in 
[23j where it was shown that w decreases monotonically 
with a, a behavior qualitatively similar to that obtained 
for w{a) in j20l and pT] for a quintessence field in a nearly 
flat potential and a scalar field rolling near a local maxi- 
mum of the potential, respectively. For < < 1, the 
behavior of w(a) also consists of a slow and monotonic 
increase in w. 

In the case of < 0, the behavior of w{a) is os- 
cillatory, and can be analytically approximated through 
a combination of sinusoidal functions [52]. For this to 
happen, we need to have T^"(0^)/F(</)^) > 3/4, requiring 
the potential to be significantly curved at its minimum. 
Here, we won't be using the analytical approximation 
for w(a) derived in [22] . but instead calculate w directly 
from equation Q by numerically solving equation ^ as 
a function of a. This will depend on three parameters: 
ri^Oj and the initial value for the field, (pi. The field's 



FIG. 1. Evolution of w for Model I, assuming different values 
ofK^: K'^ = (dot-dashed red line); K'^ = -10 (long dashed 
black line); = —20 (short dashed red line), and = —30 
(solid black line). The values of fi^^ and were fixed at 0.74 
and 0.30, respectively. 



initial velocity, (p is assumed to be zero, which determines 
the initial value of the equation of state to be —1. Also, 
we will be assuming a quadratic potential, such that 



V{<P) = p^„ + V2<P\ 



(12) 



where V2 is the curvature at the minimum, determined 
by V2 = (3/8) (1 - K"-) p^o- 

In Fig. [l]we show examples of the evolution of w{a) 
for different values of K^, with fl^Q and (pi kept fixed, 
and Hq = 70 kms^^ Mpc"^ We can see that w evolves 
monotonically for — 0. As the value of becomes 
negative and increases in absolute value, the potential 
becomes more curved at the minimum, and w exhibits 
a more dynamical behavior. For these negative val- 
ues, the field is actually able to significantly overcome the 
Hubble friction term and cross the potential's minimum, 
oscillating around it with decreasing amplitude. This is 
reflected in the equation of state in the form of maxima 
with decreasing amplitude located between points where 
the field's motion comes to a halt, which happens when 
the field reaches the maximum height it possibly can on 
the potential. As the curvature at the minimum of the 
potential, and thus K^, gets larger, the oscillatory be- 
havior becomes more evident, due to a decrease in the 
amount of time required to complete a full oscillation 
about the minimum. However, note that if (pi — 0, then, 
whichever values the other parameters might take, no 
oscillations occur, and the field effectively behaves as a 
cosmological constant. 



B. Model II 

The second model we consider is that of a scalar field 
oscillating around zero, with a potential of the form [27] 
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(13) 



This potential is very flat for small field values, acquir- 
ing a quadratic form for large field values, so that 



{(j)/M) 

(mV2) 02 



2a 



» M, 
« M. 



(14) 



The scale m determines the curvature of the potential 
at the minimum, V"{0), whereas the scale M deter- 
mines where the potential changes shape and enters the 
quadratic region, where the field will eventually oscillate. 
This model allows, if certain conditions are met, to have 
a field slowly-rolling until recently and then present some 
oscillatory behavior close to the present. 

It has been show that for a power-law potential, 
V{4>) oc |(^|", one obtains w = {n — 2){n + 2) when aver- 
aged over the oscillation period, T, as long as it is much 
smaller that the Hubble time, i.e., T « H'^ piH^ . 

Therefore, in the oscillatory region {(p << M) of the 
potential being considered, w should average to zero, be- 
cause we then effectively have V{(j)) oc i.e. n = 2. 
This means the scalar field would then be behaving like 
pressureless matter, which goes against the need for hav- 
ing a negative pressure component to drive the cosmic 
acceleration. However, this model can be fine-tuned so 
that the oscillatory region of the potential is reached close 
enough to the present that the predicted large-scale dy- 
namic of the Universe does not deviate much from that 
expected in a model where it is the presence of a cosmo- 
logical constant that makes the Universe accelerate [?7] . 
This fine-tuning can be achieved if the parameters of the 
model being considered take values around fi^g = 0.75, 
m/Ho — 1130.6, with Hq being the present-day value of 
the Hubble parameter, M = 0.002 and (l>i/M = 23.7, 
where (jji is the initial value of the field. In this particu- 
lar realization of Model II, the oscillatory behavior only 
starts around a — 0.8, avoiding the appearance of the 
gravitational instabilities that 28J have shown to be a 
problem for early rapidly oscillating models with a neg- 
ative averaged equation of state. 

In Fig. [2] we show two examples of the evolution 
of w in Model II, for M = 0.0018 and M = 0.0022, 
with Hq = 70 kms^^Mpc^^ and all other parameters 
fixed by the ratios mentioned in the previous paragraph. 
For larger M, one would expect the field to enter the 
quadratic region of the potential sooner, meaning oscilla- 
tions should start at a lower value of a. However, because 
the ratio (j)i/M is kept fixed, an increase in M will lead 
to a proportionally larger which by itself would lead 
to a delay in the start of the oscillatory period. In fact, 
it is the change in (jji which has the largest impact on the 
value of a at the beginning of such period, given that, as 
Fig. [2] shows, increasing M leads to a lower value of a at 
the start of the oscillatory period. Before entering this 



FIG. 2. Evolution of w for Model II, assuming different 
values of M: M = 0.0018 (dashed red line), and M = 0.0022 
(solid black line). The value of fi^g was assumed to be 0.75. 



period, the field's evolution is quite similar in both cases, 
slowly rolling down the potential with w almost indis- 
tinguishable from —1. As the field enters the quadratic 
region, it starts to oscillate ever more rapidly around the 
potential's zero minimum. We will then have w varying 
between ±1, with a constant averaged value of 0. Thus, 
the field in effect starts behaving like pressureless matter, 
with its energy density decreasing as a function of a~^. 
And, since V{(f>) = when the field reaches its maxi- 
mum height on the potential during an oscillation, the 
field's oscillation amplitude should decrease with a~^^'^. 



III. EUCLID 

The Euclid mission is a medium-class space mission 
expected to be launched by the European Space Agency 
in 2020. Its main objective is to help better constrain 
the large-scale geometry of the Universe and nature of 
the dark energy and dark matter components of the Uni- 
verse. More specifically, Euclid aims to enable the mea- 
surement of the dark energy parameters wq and Wa with 
2% and 10% accuracy, respectively [^S]. These parame- 
ters characterize the dark energy equation through [5D] 



w{a) = Wo + {1 — a)wa 



(15) 



where wq is the present-day value of w. 

The specifications of Euclid were defined so that its 
main scientific objective can be achieved using the infor- 
mation contained in the Baryonic Acoustic Oscillations 
(BAO) present in the matter power spectrum and in the 
Weak Lensing (WL) features of the large-scale distribu- 
tion of matter [2^. The two cosmological probes will 
rely mostly on Euclid's wide survey, although they will 
require substantially different sets of data: spectroscopic 
for BAO and photometric for WL. The Euclid wide sur- 
vey will cover an area of at least 15000 deg^, restricted to 
galactic latitudes |6| > 30°, and will be able to measure 
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the shapes and redshifts of galaxies up to a redshift of 
z « 2, with a median redshift of approximately 1 [29^ . 

A. Baryonic Acoustic Oscillations 

In the early stages of our Universe, before the re- 
combination era, its composition consisted of a dense and 
hot plasma where photons and baryons were coupled. 
This primordial plasma struggled between the possibility 
of a gravitational collapse and the repulsion of the out- 
ward electromagnetic and kinetic pressures. These two 
competing effects led to the formation of acoustic waves 
in the plioton-baryon fluid, whose signature is still im- 
printed both in the cosmic microwave background radia- 
tion (CMB) and in the large-scale distribution of matter 
(BAO). Namely, a series of peaks in the CMB angular 
power spectrum and in the galaxy power spectrum, with 
a characteristic scale, the sound horizon at recombina- 
tion, s, which can be accurately measured using present 
CMB data, such as that acquired by WMAP f^;. Its 
identification, both in the transverse, y, and radial, j/', 
directions using the galaxy power spectrum, then allows 
the comoving distance, £^nd the Hubble parameter, 

H{z), as a function of redshift, to be constrained, given 
that Ell 



represent the relative errors associated to the measure- 
ments of the transverse and radial scales, respectively. 
We note that iCosmo uses the analytical approximation 
of |36j to evaluate the expected accuracy of an experi- 
ment's measurement of the BAO scales. 

In our calculation of the constraints imposed on the 
free parameters of the two quintessence models being con- 
sidered, using the BAO method, we assumed the galaxy 
redshift distribution of Euclid's spectroscopic wide sur- 
vey to be that in [37], for a limiting flux of 3 x 10~^ 
erg s~"'^cm~^. We also assumed a median redshift of 1.1, 
a redshift measurement error of Az = 0.001(1 -I- z) and 
a redshift range of [0.4, 2.2] distributed over 14 redshift 
bins, following |29j . 



B. Weak Lensing 

Weak Lensing is a subtle gravitational lensing effect, 
caused by the slight deflection of the light emitted by 
distant galaxies due to the existing matter distribution 
along their line of sight. This will lead to changes in 
the perceived galaxy shapes. For instance, a galaxy that 
is intrinsically circular could appear to us as elliptical. 
Such type of distortion can be approximated to first order 
through a simple matrix transformation [311 El] 



and 



y 



H{z)s 



(16) 



(17) 



In order to assess the constraining power of the BAO 
method applied to Euclid's spectroscopic survey with re- 
spect to the parameters that will be varied in the two 
quintessence models we are considering, we will use the 
iCosmo software package [33] to determine the Fisher ma- 
trix associated to this survey. The Fisher matrix is a 
useful tool for planning future observations, as the diag- 
onal elements of its inverse indicate the best possible con- 
straints that can be achieved for the assumed free param- 
eters, according to the Cramer-Rao theorem [35 . The 
Fisher matrix associated with the BAO method, when 
the available information allows for the estimation of the 
two characteristic acoustic scales y and y' , is given by 
[35] 



pBAO 



dv(zi) dyjzj) 



dy'jzj) dy'(zi) 



(18) 



dp" 



dpf 



where the sums run over the observational bins at differ- 
ent redshifts z, p°' represents the parameter with respect 
to which the partial derivative is taken, while x and x' 



X2 

y2 



-K - 71 

-72 



-72 

-K + 71 



Xi 

Vi 



(19) 



The components 71 and 72 form the complex shear, equal 
to 71-1-172. A positive 71 stretches the image along the x 
direction and compresses it on the y direction, 72 has the 
same effect, but respectively along the y = x and y = —x 
directions, while k is the convergence, measuring the di- 
latation or contraction of the observed image. The pair 
of coordinates (xi^yi) and (0:2, 2/2) denote points on the 
observed image and on the source, respectively. The mag- 
nitude of the lensing signal depends both on the amount 
of matter along the line of sight, and on the distances be- 
tween the observer, the lens and the source. This makes 
weak lensing ideal for measuring the Universe's mass dis- 
tribution and geometry and, therefore, for constraining 
cosmological parameters, such as those that can be used 
to characterize the dark energy. 

We again use the iCosmo software [33] to forecast the 
constraining power of Euclid's WL survey with respect to 
the parameters that will be varied in the two quintessence 
models under consideration. Euclid will use a weak lens- 
ing technique known as tomography: the galaxies in Eu- 
clid's wide survey will be divided into separate (photo- 
metric) redshift bins, with the two-point correlation func- 
tion of the shear field computed by looking at the shear 
patterns in each bin and between different bins. The 
corresponding cross-power spectrum between bins i and 
j characterized through the coefficients Cij{l), where I is 
the order of the multipole moment. The Fisher matrix 
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associated with weak lensing tomography is then given 

by [Ml SQl E] 



I 



21 



1 dOjil} , dCkmil) 

-J sky ['■-^lljk 



(20) 



where again represents the parameter with respect to 
which the partial derivative is calculated, and fsky is the 
fraction of the sky covered by the survey. The covariance 
matrix, [C;],^^, for a given I for the i — j bin pair, can be 
written as [JT] 



'~1int r 



(21) 



where 7i„t is the random mean-square intrinsic shear in 
one component, 8ij is the Kronecker symbol and Ng is 
the number density of galaxies in the i-th bin. 

The coefficients that characterize the shear cross-power 
spectrum between redshift bins i and j, Cij{l), are given 

by 



a, {I) = 



x{z?H{z) 



X{z) 



(22) 



where only the Fourier modes with / ~ kx are assumed 
to contribute to the integral. W{z) are window functions 
(also known as lensing efficiency functions), dependent on 
geometric factors, and Pg is the three-dimensional matter 
density power spectrum at redshift z, which depends on 
how fast is the growth of the density perturbations in the 
matter distribution. 

Assuming that the Universe contains only photons, 
baryons, cold dark matter and a quintessence scalar field, 
as long as this last component does not dominate the to- 
tal energy density of the Universe at early times, namely 
before recombination, the Eisenstein & Hu transfer func- 
tion |43j can be used to determine the power spectrum of 
the initial matter density perturbations. In order to in- 
fer the non-linear corrections to the evolving linear power 
spectrum, we used the HALOFIT model [44], although 
this may lead to some loss of precision in our results 
since, as showed in [IS], use of the HALOFIT model in 
dark energy models with evolving w requires suitable cor- 
rections if one is interested in an accurate estimate of the 
non-linear corrections to the linear power spectrum. 

Following [5^ , we assume the Euclid weak lensing pho- 
tometric survey will measure the shape of over 2 billion 
galaxies at a median redshift of 0.9, covering a redshift 
range between < z < 2 distributed over 10 redshift bins 
in a 15.000 deg^ sky. For calculation purposes, the as- 
sumed number of galaxies per arcmin^ is 30. The galaxy 
redshift distribution is expected to be close to that given 
by the analytical formula of [46], with parameters a = 2 
and /3 = 1.5 [IHIHI]- The photometric redshift measure- 
ment error is assumed to be Az/(1 + z) = 0.05, while 



lint was taken to be 0.22 [IT]. The multipole range used 
in the analysis is limited to ^max = 5000 to minimize the 
effects of baryonic feedback on the power spectrum [37] . 

We should point out that, in our calculations, we have 
neglected the effect of intrinsic alignments. Hence our 
forecasts are optimistic in that sense. 



IV. RESULTS 

In this section, we present the results for Euclid's ex- 
pected constraints on the parameters of the quintessence 
models under consideration. Given that there aren't a 
priori theoretically preferred values for these parameters, 
we will set their fiducial values to those that maximize the 
joint posterior probability distribution of the parameters 
given the most recent data available from the Supernova 
Cosmology Project (SCP) [35]. 



A. Fiducial parameter combinations 

Supernovae type la are standard candles, since they 
reach essentially the same luminosity at their peak 
brightness. The ratio of this constant peak luminosity 
to the measured flux is proportional to the luminosity 
distance squared. If we also know the redshift at which 
a supernova occurs, we can then determine the relation 
between co-moving distance and redshift, which carries 
information on all the cosmological parameters that af- 
fect it. Namely, those that characterize the quintessence 
models we are considering. 

Assuming that, in each of the quintessence models be- 
ing considered, all possible values for the free parame- 
ters have equal prior probability, then the joint posterior 
probability distribution of the parameters for each model 
is proportional to the likelihood of the most recent SCP 
data. This is provided as 



/i(z) — m — M, 



(23) 



where z, m and AI are, respectively, the redshift, (mea- 
sured) apparent and (assumed) absolute magnitudes of 
the observed supernovae. The quantity fi{z) is related 
to the cosmology-dependent luminosity distance, (in 
parsecs), through 



//(z) = Slogigdz 



(24) 



Given that the uncertainty associated with the SCP 
estimates of fi{z) is assumed to be well characterized 
through a Gaussian probability distribution, the likeli- 
hood of the most recent SCP data is then proportional 
to the exponential of — x^/2, where [35] 

X' = T.iMp) - ^^^iz^)] [E-']^^ IMp) - A^,(z,)]. (25) 
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with the fii{zi) and representing, respectively, the 
SCP data and the inverse of its covariance matrix, which 
includes systematic uncertainties. The fi{zi,p) are the 
theoretically expected values for the fii, assuming the 
measured redshift Zi and the set of parameters p for the 
model / to be true. 

For the Model I, we have allowed ^^o, (f>i and 
to vary freely within the intervals [0.70,0.80], [0.0,1.0] 
and [—50.0,0.0], respectively. In the case of Model II, 
we have allowed fi^o and M to vary freely within the 
intervals [0.60,0.80] and [0.0180,0.00220], respectively. 
The interval allowed for M, as well as the conditions 

= 23.7 M and m — 1130.6 iJ^j"^, were set to force the 
scalar field to only starts oscillating around a — 0.8, ren- 
dering the model viable as discussed before. The present- 
day value of the Hubble parameter was always assumed 
to be Hq — 70 kms^^ Mpc^^. For each model, the com- 
bination of parameters that maximizes the likelihood of 
the most recent SCP data, and hence the joint posterior 
given our assumptions of flat priors for the parameters, 
is 



• Model I: n 



4>o 



0.74, K"^ = -15.0, 



0.30 



• Model II: n^o = 0.75, M = 0.00202 



which we will henceforth take to be our fiducial model 
parameter sets. We will also always assume to be work- 
ing on a spatially-fiat Universe, thus flmo = 1 — ^4>o^ 
where the power spectrum of primordial adiabatic Gaus- 
sian matter density perturbations is approximately scale- 
invariant Ug — 1, the abundance of baryons is taken to 
be 0.045 and the present-day amplitude of the matter 
density perturbations smoothed on a sphere with a ra- 
dius of 8 h~^ Mpc is equal to erg = 0.8, motivated by the 
WMAP-7 results [13 EI]. Fig. |3] shows that the recent 
evolution of the Hubble parameter, H, as a function of 
redshift z, for the models we have been considering and 
assuming the fiducial parameters, does not deviate much 
from the KCDM case. 



B. Confidence Regions for Euclid 

Using the results of last section, we forecast the con- 
straining power that the Euclid mission will have on the 
models analyzed in this work. For that, we have used 
the software iCosmo [HS], using as input the fiducial val- 
ues defined before and leaving all other parameters with 
their default values. We then marginalized over the non 
interesting parameters and obtained the two dimensional 
68% and 95% confidence regions, which we plot here, 
alongside with the one dimensional distribution for each 
parameter of interest. 

The first figure we present here, Fig. |4j shows the ex- 
pected 68% and 95% confidence regions obtained for the 
two quintessence models analyzed in this work, when the 
Weak Lensing method is applied to Euclid's photomet- 
ric survey. It is clear that through this survey it will 




FIG. 3. Recent evolution of the ratio between the Hubble 
Parameter, H{a), as a function of the scale factor a, evalu- 
ated for the fiducial parameters of both models and the corre- 
sponding AC DM (with Oa = fi^o). The black (full) and red 
(dashed) lines correspond to Models I and II, respectively. 



be possible to place tight constraints on these particular 
quintessence models. 

On the other hand, Fig. [5] shows the forecasted confi- 
dence regions for both models when the BAO method is 
applied to Euclid's spectroscopic survey. It is clear that 
through this survey it will not be possible to place as tight 
constraints on these particular dark energy models as by 
applying the Weak Lensing method to Euclid's photo- 
metric survey, given that the confidence regions are now 
much wider around the central values of the parameters. 

In the case of Model I, a clear negative correlation is 
found between the preferred values for and 0i, using 
both the WL and BAO data. This is expected given 
that a more negative would lead to earlier oscillations 
in the quintessence field, which could be prevented by 
choosing higher values for 0^. However, interestingly, the 
relationship between Jl^o and (as well as (/f^) changes 
whether its the WL or the BAO data that it is being 
considered: a positive (negative for correlation for 
the first type of data becomes negative (positive for 0^) 
for the second type of data. 

In the case of Model II, there is a positive correla- 
tion between the two parameters. In general, the oscil- 
lations in the field start sooner both for larger M and 
rj^o- However, fixing the ratio (t>i/M, to keep the start 
of the oscillatory period roughly constant, produces the 
opposite behaviour as discussed earlier on. Hence, an in- 
crease of AI (and the corresponding decrease of (pi) are 
compensated by a larger fi^o producing thus a positive 
correlation. 

Lastly, Fig. [6j shows the forecasted confidence regions 
for both models arising from the combined application of 
the Weak Lensing and BAO methods to the Euclid wide 
survey. As expected, the combination of the constraints 
obtained through the two methods considered does not 
show a great improvement over those resulting from ap- 
plying just the Weak Lensing method. Tables [l] and [H] 
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FIG. 4. It is shown the 68% (red line) and 95% (black line) 
forecasted confidence regions obtained when the Weak Lens- 
ing method is applied to Euclid's photometric survey. The 
upper panel refers to model I, and the lower panel to model 
II. For both models, we have assumed as central values the 
fiducial values determined in the previous section. We also 
plot the one-dimensional distributions for each parameter. 





Weak Lensing 


BAO 


Combined 




0.000337 


0.0312 


0.000332 




3.281 


51.734 


3.190 


4>i 


0.0267 


0.438 


0.0262 



TABLE I. The la marginalized values for the first model 
parameters obtained by applying the Weak Lensing and BAO 
methods, individually and jointly, to Euclid's wide survey 




O.eB (1.70 0.75 O.BO -1.0 0.0 _ 1.0 -100 100 



-0.003 




0.65 0.75 0.65 -0,002 0.002 0.006 

11,0 " 



FIG. 5. It is shown the 68% (red line) and 95% (black line) 
forecasted confidence regions obtained when the BAO method 
is applied to Euclid's spectroscopic survey. The upper panel 
refers to model I, and the lower panel to model II. For both 
models, we have assumed as central values the fiducial values 
determined in the previous section. We also plot the one- 
dimensional distributions for each parameter. 





Weak Lensing 


BAO 


Combined 




0.000380 


0.0350 


0.000358 


M 


0.000253 


0.00155 


0.000236 



TABLE II. The Icr marginalized values for the second model 
parameters obtained by applying the Weak Lensing and BAO 
methods, individually and jointly, to Euclid's wide survey 



summarize all the results by presenting the Icr values for 
the first and second models parameters, respectively. 

V. CONCLUSION 

In this paper we considered two quintessence models 
and determined the expected constraints that the future 
space mission Euclid will be able to place on them. The 
chosen models have the particular feature of predicting 



recently oscillating equations of state, in opposition to 
the current paradigm of AC DAI, for which the equation 
of state takes a constant value of —1 throughout cosmic 
evolution. 

In order to obtain plausible fiducial values for the pa- 
rameters of the models, we have used the most recent 
data available from the Supernova Cosmology Project. 
We performed a statistical analysis, whereby the param- 
eters were allowed to vary within plausible a priori inter- 
vals, so that the combination that maximized the joint 
posterior of the parameters, given the Supernovae data, 
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0.7395 0.7400 0.7105 0.35 0.30 D.35 -30 -15 -10 



0,00S5 
SQ.OQSO 
0,0015 




0.740 0.750 0.751 0.0015 0.0030 0.0025 



could be estimated. As expected, in both cases we ob- 
tained fiducial values for the parameters that lead to an 
overall expansion history for the Universe quite similar 
to that of ACDM, although the recent evolution of the 
equation of state in both models is very different with 
respect to the ACDM case. 

Then, using the software iCosmo, we have obtained 
the expected constraints one will be able to place on 
the quintessence models considered, around the assumed 
fiducial values, using Euclid data. We concluded that, 
while applying the Weak Lensing method to Euclid's 
photometric survey it will be possible to place tight con- 
straints on these models, the BAO method applied to Eu- 
clid's spectroscopic survey per se will not allow that. In 
order for the latter to approach the constraining power of 
the former, one will have to consider the full information 
contained in the spectroscopically-derived matter power 
spectrum, P{k), including its full shape, amplitude and 
redshift distortions, as suggested by forecasts for Euclid 
using the wq — Wa parameters for dark energy [SD] . Our 
results agree with [FT], where several oscillating dark en- 
ergy models were considered, and the conclusion reached 
was that those oscillations would probably be more eas- 
ily detectable using Weak Lensing data, rendering this 
as, possibly, the best approach to constrain the type of 
models we have been considering. 
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